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ABSTRACT 

We have used a multiwavelength data set from the Canadian Galactic Plane Survey 
(CGPS) to study the Galactic Hn region KR 140, both on the scale of the nebula itself 
and in the context of the star forming activity in the nearby W3/W4/W5 complex of 
molecular clouds and Hn regions. From both radio and infrared data we have found a 
covering factor of about 0.5 for KR 140 and we interpret the nebula as a bowl-shaped 
region viewed close to face on. Extinction measurements place the region on the near 
side of its parent molecular cloud. The nebula is kept ionized by one 08.5 V(e) star, 
VES 735, which is less than a few million years old. CO data show that VES 735 
has disrupted much of the original molecular cloud for which the estimated mass and 
density are about 5000 M and 100 cm~ 3 , respectively. KR 140 is isolated from the 
nearest star forming activity, in W3. Our data suggest that KR 140 is an example of 
spontaneous (i.e., non-triggered) formation of, unusually, a high mass star. 



Subject headings: HII regions — ISM: individual (KR 140) — stars: formation 



1. INTRODUCTION 

Massive OB stars are almost always found in clusters. In fact, accumulating observational 
evidence suggests that most stars, regardless of mass, actually form as members of some kind of 
group, cluster, or association. While there has been a lot of work to understand the processes 
involved in the formation of a single star (e.g., Shu, Adams, &; Lizano 1987), theories of cluster 
formation are still in their infancy (see the recent reviews by Elmegreen et al. 2000, and Clarke, 
Bonnell & Hillenbrand 2000). Nevertheless, many important issues have already been identified. 
Foremost among these is the question of whether the formation of clusters, particularly ones with 
OB stars, are triggered by an external agent (Elmegreen 1992) such as expanding Hn regions 
(Elmegreen h Lada 1977), or colliding molecular clouds (Loren 1976, 1977; Scoville, Sanders & 
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Clemens 1986; Usami, Hanawa, & Fujimoto 1995). Despite triggered or sequential star formation 
being theoretically and intuitively appealing, a major problem is to determine unambiguously 
a cause and effect relationship, because of the long time scales (and time lags) of the processes 
involved. There is observational evidence for regions that have had triggered star formation, both 
on large scales (e.g., IC 1396, Patel et al. 1998) and on small scales (e.g., IC 1805, Heyer et al. 
1996). However, there are other young star forming regions, like Taurus, where no evidence of a 
trigger can been found. These are often in modest-sized molecular clouds and contain only lower 
mass stars. 

The Perseus Arm star forming regions W3/W4/W5 (Westerhout 1958) have been studied 
extensively over the last twenty years (e.g., Lada et al. 1978; Braunsfurth 1983; Digel et al. 1996; 
Normandeau, Taylor & Dewdney 1997; Heyer & Terebey 1998). They are often considered to be 
the archetypical examples of how the formation of massive star clusters can be triggered by the 
influence of other nearby clusters. For example, W3 is thought to be have been triggered by the 
expansion of W4 (Dickel 1980; Thronson, Lada, & Hewagama 1985; van der Werf & Goss 1990) 
and there is evidence that the expansion of W5 is also triggering star formation (Vallee, Hughes, 
& Viner 1979; Wilking et al. 1984). 

In new high resolution multiwavelength (radio and mid-infrared) data of the W3/W4/W5 
complex from the Canadian Galactic Plane Survey (CGPS; English et al. 1998) we have identified 
a star forming region containing a single O star. Figure Q shows a 1420 MHz continuum image 
from the CGPS pilot project (Normandeau et al. 1997). The circled area is the Hn region in 
question, KR 140 (I = 133.425°, b = 0.055°; Kallas & Reich 1980). This region appears to be 
completely separate from the vigorous star formation going on nearby in W3, although it is in the 
same Perseus arm molecular complex (see context in CO images in Heyer et al. 1998). What is 
unusual is that this massive star seems to have been formed spontaneously. 

In this paper we present and analyze the multiwavelength data on KR 140 in order to quantify 
the properties of this region of spontaneous massive star formation. As described in § [2] these 
data have sufficient resolution (1') to resolve this Hn region for the first time and reveal a fairly 
symmetrical structure. Based on first impressions, we thought that KR 140 might prove to be a 
"textbook" spherical Hn region. Instead, we find that KR 140 is likely a bowl-shaped Hn region 
(an example of a blister geometry: Israel 1978, Yorke et al. 1989; § || and § 
found that the Hn region is kept ionized by an 08.5V(e) star, VES 735, and is at a distance of 
2.3±0.3 kpc from the Sun (Kerton et al. 1999). 

We analyze the ionized component of KR 140 in § || and make various estimates of the age 
in § |E[ The dust, molecular, and atomic components of KR 140 are examined in § ||, § |^, and § [8|, 
respectively. In § [|| we discuss star formation in KR 140 in the context of the Perseus arm and the 
possible accompanying cluster. 



4.2.2, 



We have 
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2. THE MULTIWAVELENGTH DATA 

KR 140 was initially cataloged in a 1420 MHz radio continuum survey of the northern 
Galactic plane with the Effelsberg 100 m telescope by Kallas & Reich (1980). At the available 
resolution of 9' KR 140 was barely resolved, with a reported diameter of about 11'. Although 
KR 140 was measured in other subsequent single-dish surveys (§ |4.2| ) it had never been examined 
with a radio interferometer to provide the angular resolution necessary for the present study. 

The data analyzed here include 1420 MHz (A 21 cm) and 408 MHz (A 74 cm) continuum 
images from the Dominion Radio Astrophysical Observatory (DRAO) Synthesis Telescope (Roger 
et al. 1973; Veidt et al. 1985). An Hi 21-cm line data cube is also available, with spectral resolution 
2.64 km s _1 and channel spacing 1.65 km s _1 . These DRAO data were obtained during the 
pilot project of the CGPS (Normandeau et al. 1997), and made into a full 8° x 5° mosaic of the 
W3/W4/W5 star forming regions. 

We also made use of a CO (J = 1 — ► 0, 115 GHz, A 2.60 mm, spectral resolution 0.98 
km s _1 ) data cube from the Five Colleges Radio Astronomy Observatory (FCRAO) obtained in 
the complementary survey described by Heyer et al. (1998). To examine the dust components we 
processed Infrared Astronomical Satellite (IRAS) data to make HIRES mosaics at 12, 25, 60, and 
100 /im. Independent processing of the entire Galactic plane data at 60 and 100 /im has been 
released as the IRAS Galaxy Atlas (IGA, Cao et al. 1997) and we have completed a complementary 
project at 12 and 25 fim (Mid-Infrared Galaxy Atlas, MIGA, Kerton & Martin 2000). 

An important feature for subsequent analysis with these data sets is the common relatively 
high angular resolution achieved. The DRAO synthesis telescope images have a resolution 
of 1' x 1.14' at 1420 MHz (the 408 MHz DRAO data have proportionately 3.5 times lower 
resolution), and the CO images are at a resolution of 50" (beam sampled). The HIRES images 
have non-circular beams of size about 2' at 100 /im, 1' at 60 /im, and somewhat less at 12 and 
25 /im. For detailed intercomparison, pairs or groups of images have been convolved to the same 
beam shape. 

3. A MODEL FOR THE 3-DIMENSIONAL MORPHOLOGY 

Analysis of an Hn region benefits from a knowledge of its three-dimensional structure, but 
often the observed two-dimensional morphology is too complex to interpret. Fortunately, this is 
not the case for KR 140. Figure || shows the four HIRES images of KR 140 at 12, 25, 60, and 
100 /im. As discussed in § ||, the 100 /im image best describes the dust distribution in KR 140. 
The 100 /im HIRES image (Figure Qd) shows that most of the dust in the nebula has been swept 
into a shell-like structure and that we are observing a region which is close to being circularly 
(axially) symmetric. In the 1420 MHz image (Figure ||) KR 140 is also seen to be fairly circular 
with a central depression. As highlighted by the contour map, there are three high intensity areas 
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- two "eyes" and a "mouth" . However, where the "nose" would be is a lower intensity area. The 
simplest interpretation appears to be a three-dimensional structure with a central hole (§ |5.2.2| ). 
We have identified VES 735 as the exciting star of this Hn region (Kerton et al. 1999); it lies 
projected close to this lower intensity area of Hn emission, though not coincident with it. 

Since we have also classified the exciting star, we can draw further conclusions about the 
geometrical structure of KR 140 from a global analysis (§ 5 in Kerton et al. 1999). Even after 
taking into account the emission from dust that was not observed by IRAS, we find that the 
total luminosity of warm dust in the nebula is far less then the bolometric luminosity of the 



ionizing star VES 735 and that the covering factor of the dust is only about 0.4 - 0.5 (see § 6.4). 
Furthermore, the total radio flux at 1420 MHz is lower than expected for an ionization bounded 
nebula surrounding VES 735, with a similar implied covering factor (§ |3~2~1 ). This together 



with geometrical information to be discussed (e.g., § 4.1) implies that KR 140 is a bowl-shaped 
Hn region (cf. Roger & Irwin 1982), rather than a classical Stromgren sphere, a not unexpected 
morphology when the O star is close to the edge of its parent molecular cloud (Yorke et al. 
1983). This model is perhaps the simplest geometry that is consistent with the derived covering 
factors, but our data cannot rule out other, more complicated geometries, such as a broken shell 
or shredded and dissipated Stromgren sphere. Observations of a champagne flow (Tenorio-Tagle 
1979), perhaps via Fabry-Perot imaging, could help establish that KR 140 is truly a blister region. 

The circular symmetry observed in the dust shell (and in the ionized gas) implies that we 
are observing KR 140 almost face on (i.e., the opening of the bowl is oriented almost along the 
line-of-sight). For example, see the simulated radio maps of the R2 model by Yorke et al. (1983). 
The CO data and extinction measurements to VES 735 both suggest that the molecular cloud is 
behind the star, with the inferred opening toward us (§ |7.1[). 



4. KR 140 FROM THE PERSPECTIVE OF THE RADIO CONTINUUM 

4.1. Size and Structure 

The basic radio morphology appears to be that of a limb-brightened hemispherical shell (most 
of the free- free emission would be from the dense and thin ionization-bounded zone). The last 
contour shown (at Tf, = 5.0 K) in Figure || has a diameter of 8.5' which for a distance d = 2.3 kpc 
corresponds to a physical diameter of about 5.7 pc. It is expected that neutral material could 
extend well beyond the Hn region, and indeed there is both dust (§ |6.2|) and CO gas (§ |7.1|) around 
KR 140. The free-free surface brightness falls off most quickly (the radio contours are more tightly 
spaced) in the west-northwest part of the nebula. A sharper ionization front would occur in more 
dense (pre-existing) material. 

Other internal detail might be interpreted as variations in distance between the star and the 
background ionization front, with lower surface brightness features originating from sectors further 
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from the star; such an interpretation of radio surface brightness has been used to construct a 
topological map of the ionized surface of the Orion molecular cloud behind the Trapezium (Wen & 
O'Dell 1995). If we use the 1420 MHz surface brightness at the pixels near the projected position 
of VES 735 and the known properties of this 08.5V exciting star, the distance that the star must 
be from the back of the nebula turns out to be only 2.14 pc. Together with the above estimate of 
the diameter across the line of sight, this is consistent with the idea that KR 140 is similar to a 
hollowed-out hemispherical bowl. 



4.2. Physical Properties 

In this section we use the 1420 MHz data to derive a number of physical properties of the 
nebula. For ease of reference, Table |I| summarizes these properties, some of the properties of the 
exciting star, and other quantities that we later derive from the infrared and CO data. 



4-2.1. Emission Measure, Radio Flux, Electron Density, and Ionized Mass 

The KR 140 emission is optically thin at 1420 MHz with a peak brightness temperature about 
10 K. For an optically-thin nebula, the specific intensity (surface brightness) I u and brightness 
temperature for thermal radiation are (Osterbrock 1989): 

I u = 2v 2 kT bu /c 2 = j u J nin e ds = j u E, (1) 

where j u is the free-free emissivity, ni is the ion density, n e is the electron density, ds measures 
distance along the line of sight, and the integral E is called the emission measure. For j v of a pure 
hydrogen nebula (a good approximation; see § |4.2.2 for models which include ISM abundances), 



E = 5.77 x 10 2 T 6! ,(T e /7500K) a35 (z7l420MHz) 2 - 1 cm' 6 pc. The peak emission measure in the map 
is about 6000 cm -6 pc and after background subtraction (E) = 2050 cm~ 6 pc. After subtracting 
the local Galactic background we find the total flux {F u = f I u dVt) to be 2.35 ± 0.05 Jy at 
1420 MHz. This is in good agreement with the value 2.3 ± 0.2 Jy reported by Kallas & Reich 
(1980). Becker, White, & Edwards (1991) detected this Hn region in a 6 cm survey with the 
NRAO 91 m. Their source, designated [BWE91]0216+6053, had F v = 2.5 Jy. Using the same 
telescope and frequency Taylor & Gregory (1983) and Gregory & Taylor (1986) record this source 
as GT 0216+608 in their survey; in their analysis they treated all sources as point sources and so 
find systematically low fluxes for extended nebulae like this (648 mJy for KR 140). 

In the optically thin limit the theoretical expectation for the spectral index for thermal 
radiation (defined as I u oc u a ) is a pa —0.1 (Oster 1961; Gordon 1988). The spectral index was 
measured between 1420 and 408 MHz. First, the distribution across the nebula was determined 
using L. Higgs' "specmap" routine (Zhang & Higgs 1997), which evaluates the spectral index for 
the spatially variable component of the emission (this approach has the advantage of circumventing 
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uncertainties in background subtraction in each image); for this procedure, we convolved the 
1420 MHz image to the lower resolution of the 408 MHz image. The typical value of a is 
0.02 ± 0.08. We confirmed this by measuring the total background-subtracted flux at 408 MHz 
as well. Considering there is a 10% uncertainty in the flux calibration in the 408 MHz pilot 
project data at the position of KR 140 (Taylor 1999), the derived spectral index is in satisfactory 
agreement with theory. 

Depending on details of the geometry, the appropriate pathlength for estimating the density 
would be of order the observed radius; with this choice we find an rms n e = 27 cm~ 3 . The mass 
of radio-observable ionized gas (allowing for He) is 163 (27 cm~ 3 /n e ) Mq. For comparison, the 
exciting star of KR 140 has a mass of 25 Mq (Kerton et al. 1999). 



4-2.2. Number of Ionizing Photons 

In ionization equilibrium the number of ionizations (where H is the predominant species) 
that occur each second in the nebula is equal to the number of (H + ) recombinations each second, 
both locally and globally. An ionization-bounded nebula occurs when there is sufficient material 
to intercept all of the ionizing photons (A < 912 A), Q(H°), emitted by the star per second 
(Osterbrock 1989). Formally, the spatially-integrated radio emission of an ionization-bounded 
nebula is an effective calorimeter for Q(H°): 

fQ(H°) = J n e n R+ a B dV = a B F u d 2 /j u (2) 

where as is the case B hydrogen recombination coefficient of H + (~ 3.3 x 10~ 13 cm 3 s _1 at 
7500 K; Storey & Hummer 1995), dV is the volume element, and we have used dV = dsdQ x d 2 
and Equation ([!]) to find the right hand side. Here / is the product of a number of correction 
factors including / c j, the all-important covering factor by an ionization front: if the star is only 
partially surrounded by gas then a fraction of the ionizing photons will escape and a lower F v will 
result. 

Equation |2] is based on recombination and emissivity for a pure H nebula at a fixed 
temperature and so / is made up of a number of factors which modify the expected radio 
continuum flux for a given Q(H°). We used the spectral synthesis code Cloudy (Ferland et al. 
1998) to analyze the observations of KR 140 in the radio, particularly the correction factors which 
comprise /. We constructed a number of toy models for differential comparisons, using a 35,500 K 
blackbody with log[Q(# )] = 48.06 and n e = 50 cm 3 ; see Table |2[ The differential results are 
appropriate for any low-density Hn region heated by a late O-type star. 

The temperature of the nebula enters in the equation through the as and j v factors, although 
the effect on the ratio is very small. Comparing models 1 and 5 we estimate fx = 0.99. 

The effect of the addition of other elements to the nebula is to increase the total radio flux. 
This effect is primarily due to the presence of He + increasing the effective free-free emissivity (He 
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has a tiny effect on the amount of ionization of H; e.g., Osterbrock 1989). Comparing models 1 
and 4, or models 2 and 5 we estimate that fne = 1-07. 

Dust competes with the gas for the absorption of ionizing photons and thus, when present 
in the nebula, will reduce the observed radio F v for a given Q(H°). Comparing models 1 and 2 
(with grains) to model 3 (without grains) we derive fdust = 0.78 for the typical ISM grains used in 
model 1 and fdust = 0.93 for dark-cloud dust, similar to the dust found in Orion, used in model 2. 
The dark-cloud dust is not as efficient as typical ISM dust in absorbing short-wavelength photons 
and so has less of an effect on the emergent radio flux. 

The combination of these factors, except for the covering factor, leads to a factor of 
f model = 0.91 ± 0.08 depending upon the grain composition (i.e., / = 0.91/ c j). 

In our study of the exciting star VES 735 (Kerton et al. 1999) we used the measured radio 
flux at 1420 MHz (2.35 ± 0.05 Jy) to calculate log[fQ(H )} = 48.05 ± 0.11, and then estimated 
f ci w 0.4 - 0.5 based on the known spectral type of VES 735 (08.5 V(e); \og[Q{H )] ~ 48.45, 
Panagia 1973). We have not attempted to put any formal error estimate on the covering factor, 
but given the uncertainties involved, both in the stellar properties and in estimating the various 
/ factors, this result is certainly consistent with the idea of KR 140 being an open bowl-shaped 
region. As demonstrated in § ^T4|, an analysis of the IR data leads to the same conclusion 
independently. 

5. THE AGE OF KR 140 

The approach we have adopted here is first to date the Hn region using data we have on the 
exciting star VES 735. With that age in mind, we then investigate the dynamics of the region, the 
goal being to show that certain scenarios for the evolution of KR 140, such as it being a blister 
region, are at least consistent with the age suggested from the exciting star. This approach is 
similar to that adopted by Dorland et al. (1986) in their study of the Rosette Nebula. 

5.1. Stellar Content 

The idea of using the stellar content of a Hii region to measure the age of the nebulosity was 
first attempted by Hjellming (1968). Basically one plots the evolutionary tracks of stars with 
various masses in the log(L/L ) vs. logT e ff plane. One obtains T e ff from the spectral type of 
the exciting star and log(L/ L & ) from the radio flux (much like Q(H )) and then the position in 
the plane determines an age for the star and thus the Hii region. Clearly the effectiveness of this 
technique depends strongly on the quality of the calibration between theoretical and observational 
quantities as well as the quality of stellar models, which have vastly improved in the thirty years 
since this technique was introduced. In the early work the primary result was to indicate whether 
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the Hn regions were ionization or density bounded, and that many Hn regions required additional, 
unobserved, sources of ionization. 

Here we follow a slightly different technique compared to Hjellming (1968) in order to avoid 
any uncertainties associated with the covering factor and structure of the nebula. We instead 
use the absolute magnitude (My) as a measure of log(L/L Q ), which is possible because we have 
a good estimate of the distance to the star and the extinction along the line of sight (Kerton et 
al. 1999). Figure |] plots the stellar evolution models of Schaerer & de Koter (1997) for a 20, 25 
and 40 Mq star along with the value determined for VES 735. The observed values for VES 735 
are consistent with a 25 M Q star with an age of a few million years away from the ZAMS (see 
Figure Q), where the age would be only of order 10 5 years. Of course, the spectral type alone gives 
us a simple upper limit to the age of KR 140: for an 08.5 V star the main sequence lifetime is 
~ 6 x 10 6 (Chieffi et al. 1998). 



5.2. Dynamical Models: Spherically-symmetric 

5.2.1. Stromgren Sphere 

The simplest description of a Hn region is that of the formation and expansion of a ionized ball 
of pure hydrogen at constant temperature in an uniform medium of constant density (Stromgren 
1939). We summarize this only as a point of departure and contrast. The evolution of a Stromgren 
sphere starts with a formation phase where the O star ionizes a region of space around it to the 
radius (R s ) given by: 

The initial rapid expansion to a radius of r, occurs in time t = (— l/n e a.B) ln(l — (rj/i? s ) 3 ) (valid 
to ri/R s ~ 0.98, Osterbrock 1989). The time to create the initial Stromgren sphere is about 
10 5 /n e (cm -3 ) y, instantaneous compared to other timescales. For KR 140, the observed radius 
(R obs ) is 2.85 pc. Taking log[Q(H )] = 48.45 for an 08.5 V(e) star like VES 735 (Panagia 1973) 
and nn 2 = 100 cm -3 for the original molecular cloud, R s = 1.2 pc. The initial Hn region is 
overpressured compared to the surroundings and will expand into a uniform medium according to 
a i 4 / 7 law: 



where C/j is the sound speed in the ionized medium, r is the initial radius and is the radius at 
time t. With the densities quoted above the Hn region will evolve from R s to R bs on a timescale 
of 10 5 years, which is improbably short. Increasing the initial density makes R s smaller and forces 
the pressure expansion stage to be longer. A value of n# 2 = 500 cm -3 will push the timescale to 
10 6 years. However, this value is not consistent with our observations of the molecular material 
(§ 0), so the Stromgren sphere is not an appropriate dynamical model for KR 140. 
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5.2.2. Stellar Winds 

In the radio image a local minimum is evident near the position of the central star (see 
Figure 1 in Kerton et al. 1999). One interpretation is that this is a wind blown bubble around the 
08.5 V(e) star. The apparent radius of the central hole is 1.4 pc. We used the model of Castor et 
al. (1975) for the size of a circumstellar shell: 

R = 28 |^^oo y /5 t 3/ 5pc (5) 

where R is the shell's radius (pc), Mq is the mass loss rate (10 -6 M & yr _1 ), V2000 is the wind 
velocity (2000 km s _1 ) n is the gas density (cm -3 ), and t% is the time (10 6 years). Stellar wind 
properties of the B2 model of Schaerer & de Koter (1997) were used. For a wide range of densities 
we obtain timescales on the order of only 10 4 years. This is far too low a timescale and clearly this 
model is not an appropriate description of VES 735. Either the mass loss rate adopted is much 
too high, or the assumption of spherical symmetry is not valid. 



5.3. Dynamical Models: Open Geometry 



Evidence summarized thus far (§ [4.2.2 and § |6.4j ) indicates an open geometry with a covering 



factor of about 0.5 for both dust and ionized material. In this case the high pressure in the interior 
of the initial Stromgren sphere or stellar wind bubble is relieved, slowing down the expansion. 
Models have been developed to investigate the evolution of an H11 region at the edge of a molecular 
cloud (see review of early models by Yorke (1986), and Comeron (1997) for a recent example). 
While complex numerical modeling is possible, Franco et al. (1994) have shown that the expansion 
of the ionization bounded side of such a blister Hn region is well described by 

ri/r D = 1 + — ■ (6) 



2r ) 

This is based upon mass conservation between material being ejected in the flow and molecular 
material being eroded off the cloud. One can envisage the evolution of a blister as consisting of 
three stages: the initial rapid formation stage, a pressure driven expansion stage, and finally a 
blow-out stage. The relative length of time of the latter two stages depends upon the distance of 
the star from the edge of the cloud and the density structure of the cloud. One very important 
point is that an O star very close (~ 1R S ) to the edge of a cloud will very quickly develop a 
covering factor of ~ 0.5 in order 10 5 years and will maintain this covering factor over the lifetime 
of the O star (Yorke et al. 1983, 1989). We assume that the star formed very close to the edge of 
the cloud, thus ignoring the pressure driven expansion stage. Using Equation |j the age of the 
region is of order a million years when R s = 0.64 pc. For the stellar properties of VES 735 this 
requires njj 2 ~ 280 cm -3 . This is somewhat encouraging as it does not require as vast a difference 
between the properties of the observed molecular cloud and the putative initial conditions. 



-10- 



6. KR 140 IN THE INFRARED 

6.1. Morphology and Emission Mechanisms 

Stars that are forming and evolving in the ISM interact with the interstellar dust component 
by heating, redistributing, and possibly destroying it. The dust can be heated by at least three 
distinct mechanisms: direct radiation from the central star, reprocessed radiation from the ionized 
gas, and diffuse radiation from the interstellar radiation field. Radiation pressure from the 
star acts on the dust in the ionized zone (Spitzer 1978) which causes the dust (and gas) to be 
pushed away from the star. Some forms of dust like poly cyclic aromatic hydrocarbons (PAHs) 
are destroyed in intense ultraviolet fields. A morphological study of the infrared emission from 
KR 140 is therefore important to fully understand the energetics and effects on the environment. 

As mentioned in § g IRAS scans of KR 140 at 12, 25, 60, and 100 ^m were processed by the 
HIRES software to generate maps of about V resolution. The new beam shapes are somewhat 
elliptical and so point sources will be visibly stretched; however, the larger scale morphology of 
the observed dust emission from the nebula will not be greatly affected by the asymmetric beam. 

The intensity of the dust emission, I dust ■> a t a particular frequency, u, from a distance 
increment ds along the line of sight has the form 



Ldust 



Ndust^a 2 B y {T dust )Q lJ (a, T dust )ds, (7) 



where N dus t is the number density of grains, B v (T,iust) is the Planck function at a temperature 
Tdust, a is the radius of the particles, and Q v (a, T dust ) is the absorption efficiency factor. 

Figure [2] shows the four HIRES images of KR 140, with the 12 and 25 ^m images convolved 
to the 1420 MHz resolution, and overlaid by 1420 MHz contours (Figure [2]a & b). The 12 and 



25 jxm emission is spread well outside the radio contours of KR 140 (this is taken up in § 6.3). 
Features that are common to all four images in Figure [2] are the bright arcs on either side of the 
nebula (easiest to see in Figure §c) that extend outside the radio contours. We interpret this as 
the limb-brightened warm dust shell around KR 140. Note that there could be cooler dust further 
out around the KR 140 complex that will not have been detected in the IRAS bands. We account 
for the energetics of this cooler dust in our models of KR 140 (§ 



6.2. Temperature and Column Density 

Equation |7| can be used to calculate a mean dust temperature for each pixel of IR emission. 
Taking the intensity ratio for any two frequencies, v\ and z/2, yields: 

'--"*-r(££4b (B) 
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where T is the Tdust of Eq. [?], and the (^2/^1 ) 3+ ^ factor is from the frequency-dependent part 
of Q v (a,T), where (3 depends on the type of dust. The most common components proposed 
for interstellar dust are silicate and graphite (or some related carbonaceous material), which, at 
long wavelengths, have (3 ~ 2. This value of (3 is close to what is generally observed in the ISM 
(Lagache et al. 1998). 

The 60 /im and 100 yum data were used to calculate the dust temperature map, as they are 
the bands where classical grain emission dominates. It is likely that non-equilibrium heating of 
very small dust grains (VSGs) contributes some of the observed 60 p flux and so the derived 
temperatures are probably slight overestimates of the true grain temperatures (Boulanger et al. 
1988). The images were brought to the same resolution and background subtracted and the IPAC 
analysis program 'cttm' was used to compute a dust temperature map. An azmuthially-averaged 
radial cut of this map is shown in Figure ||. The temperature distribution was sampled every 1' in 
radius and every 10° azimuthally. The temperatures in the central region are around 31 K. Nearer 
the edge of the Hn region, the dust temperatures drop to around 28 K. However, a line of sight 
passing near the star also has cooler dust in the background and the foreground, which would 
lower the apparent temperature observed. 

A calculation of dust temperature from first principles was done as a check on the empirical 
values output from 'cttm'. In this calculation, a single silicate or graphite dust grain with a radius 
of 0.1 //m (a typical interstellar size, see Kim, Martin & Hendry (1994)) was placed at a distance 
of 3.0 pc from the center of the nebula, corresponding to a dust grain within the KR 140 dust 
shell. The luminosity of the 08.5V(e) exciting star is about 10 5 L Q (Panagia 1973). Making use 
of Planck-averaged absorption factors from Laor and Draine (1993), we find that if tjjv ra 1> the 
calculated dust temperature is about 28 K, in agreement with the temperature deduced empirically. 
Note, however, that this calculation did not take into account the recombination-line photons 
emitted by the ionized gas. Furthermore, some of the free-free photons and collisionally-excited 
cooling lines are emitted not in the ultraviolet, but in the optical (Osterbrock 1989), where the 
dust absorption efficiency is somewhat lower. 

The optical depth r v from dust is defined to be 

r v = J N dust Tra 2 Q u (a, T dust )ds, (9) 

and can be calculated by dividing I u (Equation |7|) by the Planck function B V (T), assuming a 
constant T along the line of sight, adopted from the temperature map. Figure |B] displays a dust 
optical depth map at 100 //m. The values range from about 0.0004 in the middle of the nebula to 
about 0.002 in the bright northwest rim. These values show that the dust is transparent to its 
own 100 //m emission. The minimum in the center of the nebula and the ring-shaped appearance 
is most simply interpreted as limb brightening in a thick shell of dust with the highest column 
density along the northwest rim. The latter is the same region of the nebula where the 1420 MHz 
radio contours fall off most steeply (§ [4.1[) , and where there is no CO emission (§ |7. 1| ) . 
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The optical depth in the ultraviolet lyy can be gauged using extinction curves from the 
literature. We are interested in the radial as opposed to line of sight optical depth. Judging the 
thickness of the shell from one of the arcs in Figure ^ gives a path length of about 7 x 10 18 cm. 
If we assume the dust is associated with gas at a molecular density of 100 cm -3 , the molecular 
column density of this gas is about 7 xlO 20 cm -2 . Using the extinction curves of Kim, Martin & 
Hendry (1994), we find that r v ~ 2.8 at 1100 A (far in the ultraviolet), whereas at 5500 A (in the 
optical), t v ~ 0.6. The estimated optical depths show that the dust in the arcs has a high enough 
radial optical depth to absorb most of the incident ultraviolet photons. However, since KR 140 
has a covering factor of 0.4 - 0.5 (§ [4.2.2 and § |6.4j ), the total infrared luminosity of the nebula 
will be correspondingly less than the bolometric luminosity of VES 735 (§ |6.4[ ). 



6.3. The 12 and 25 fiva Emission 

It has been known for over a decade that there is excess emission within the 12 fim IRAS 
passband (e.g., Boulanger, Baud & van Albada 1985). Onaka et al. (1996) found that more than 
70% of the 12 /iin diffuse interstellar emission detected by IRAS is emitted in spectral features 
attributed to polycyclic aromatic hydrocarbons (PAHs; Leger & Puget 1984; Allamamdola, 
Tielens & Barker 1985). The 12 fim image of KR 140 is very instructive. There is very little 
12 |im flux within the radio contours of KR 140, implying that PAHs are destroyed there in the 
intense ultraviolet radiation field. However, there is a large amount of diffuse flux outside the radio 
contours, especially to the west (Figure §a). The PAH emission is a tracer of the photodissociation 
region around a nebula (Giard et al. 1994; Bregman et al. 1995; Fig. |To|) . 

Unique among the HIRES images of KR 140, the 25 //m image (Fig. ^o; Fig. 0) shows a bright 
spot near the center of the nebula, on a transition between a radio peak (the "left eye") and the 
deepest depression. The IRAS Point Source Catalog (Joint IRAS Science Working Group 1988) 
lists this feature as IRAS 02165+6053, and it also has been identified with VES 735 (Bidelman 
1988). Figure |?j shows that IRAS 02165+6053 and VES 735 are practically coincident. Dust closer 
to the star would tend to be warmer, contributing to the spot if not pushed away by radiation 
pressure. 

The entire nebula looks hotter and more extended at 25 /jm than it would be for equilibrium 
emission from normal-sized grains. Most of the 25 ^m emission is probably contaminated by 
non-equilibrium emission from very small grains (VSGs; Sellgren 1984) which have absorbed a UV 
photon and have had their temperatures instantaneously rise to ~ 10 3 K. Although these small 
grains make up a tiny fraction of the mass in the dust distribution, they make up a good fraction 
of the number distribution and absorb a significant fraction of the near-ultraviolet radiation. This 
VSG emission is therefore a major source of uncertainty when analyzing the 25 fim image J 3 ] Note 



3 Further evidence for contamination in the 25 /im IRAS passband was presented by Cox (1990). He suggests that 
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that the "temperature-spiking" phenomenon also occurs outside the ionized zone, far from the 
exciting star. 

6.4. Infrared Models 

Hn regions are some of the most luminous objects in the Galaxy when observed in the 
infrared, especially at wavelengths longer than 60 fim. In fact, if the dust shell covers 4-7T steradians 
around the exciting star, the infrared luminosity should be a good measure of the star's bolometric 
luminosity. To account for the true extent of the dust, we can define a covering factor of dust (fed)- 
Often it is simply assumed that Lj r = L^; however, this is not correct even for f cc i = 1. Some of 
the radiation from the star and nebula is at long enough wavelengths to avoid being absorbed by 
the dust. Using Cloudy we found that for models with fed = 1, log (L^/L^i) ~ —0.1. 

To calculate the integrated infrared flux from KR 140, the background in each of the four 
HIRES images was fitted and subtracted, and the total flux from the nebula was measured in 
each band. In order to estimate L^oi using these data, we used Cloudy to simulate the emission 
from large classical grains, roughly matching the fluxes at 60 and 100 /im. This approach allows 
us to account for emission from grains with a range of temperatures and thus unobserved emission 
at long wavelengths. First, we constructed models with f C( i = 1. We find \og(Lpni) = 37.90. 
Since the resulting spectrum misses most of the observed 12 and 25 /im flux (which is caused by 
temperature spiking of VSG's and the excitation of PAH molecules) we converted the observed 12 
and 25 /im fluxes to luminosities using tophat approximations to the IRAS passbands (Emerson 
1988) and added the results to Lpm- With this addition we find \og(Lm) = 38.00. Correcting 
from Li r to we find log (L^/Lq) = 4.52. This is significantly below what would be expected 
for any late O main sequence star (e.g., Panagia 1973). We interpret this low apparent Lf, ; as 
being due to a covering factor of 0.4 - 0.5. 

We know that VES 735 is an 08.5 V(e) star. Figure || demonstrates that one can reproduce 
the observed Lj r using the appropriate stellar parameters and a covering factor of about 0.5. A 
single temperature (T = 28.25 K) v 2 B v spectrum is shown for comparison; note the wider model 
curve caused by the range of dust temperatures contributing to the emission. 

Using an alternative calibration of stellar parameters which includes wind-blanketed models 
(Kerton 1999), we find f c d ~ 0.5 and f c i ~ 0.7. These results are still consistent with a blister 
model for KR 140. 

7. KR 140 IN A MOLECULAR CLOUD 



an iron oxide emission line falls within this IRAS passband which would result in even greater excess 25 /im emission. 
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7.1. CO Signature 

A slit spectrum of VES 735 allowed us to measure the radial velocity of the nebular Ha line 
within 30" of VES 735 to be —46 ± 2.1 km s _1 with respect to the Local Standard of Rest [LSR] 
(Kerton et al. 1999). The differential radial velocity between VES 735 and the nebular line was 
measured to be +2.0 ± 2.2 km s^ 1 . 

An expected effect of the evolution of the Hn region is photodissociation of molecular gas 
both inside the Hn region and in the immediate surrounding interstellar medium (ISM). The 
CO data traces the molecular gas content of the ISM, so an examination of the CO data cube 
ought to reveal a lack of emission near the ionized gas velocity of KR 140. Indeed, a distinct CO 
hole was found within the radio contours of KR 140 over the velocity range —45.53 km s _1 to 
—47.16 km s -1 (LSR), which corresponds to three channels of width 0.813 km s -1 in the CO cube. 
At more negative velocities, CO emission from the parent molecular cloud fills in the 1420 MHz 
contours. 

Figure ^ shows the sum of these three channels overlayed with the 1420 MHz continuum 
contours. A well defined ring structure is clearly seen. Interestingly, the ring does not extend all 
the way around the nebula: there is no CO emission in the north-west. This is the same area 
where the 1420 MHz contours are falling off more sharply (§ |4.l[) , there is a bright infrared arc 
(§ |6.1[) , and a bright Hi feature (§ |8|). The relationship between 12 /im emission, which is a good 



tracer of the PDR, and the CO emission is shown in Figure |10[ The north-west peak in the 12 /im 
emission corresponds to the region where there is no CO emission. To the east, as one moves away 
from the Hn region the 12 //m peak occurs first followed by the peak in the CO emission. 



7.2. Density and Mass of the Molecular Cloud 

To estimate the mass of the molecular cloud we integrated the CO cube over the whole velocity 
range of the cloud (—45.5 km s _1 to —52.8 km s _1 ), and measured the total surface brightness of 
the molecular cloud. We can use the empirical X factor to convert from CO surface brightness to 
molecular hydrogen column density. We have used X = (1.9 ± 0.5) x 10 20 cm~ 2 (K km s -1 ) -1 
(Strong & Mattox 1996), where the quoted uncertainty is to take into account the range 
of X values measured for a variety of different clouds, to obtain a mean column density of 
N(H2)= 2.2 x 10 21 cm" 2 . Assuming a path length of 7 pc (the north-south radial extent) the 
density of the molecular cloud is roughly n# 2 = 100 cm -3 , typical of a giant molecular cloud (Blitz 
1993). 

By integrating the column density over the face of the cloud, we estimate the mass of the 
cloud to be (4.4 ± 1.6) x 10 3 M Q . This mass includes a He correction factor of 1.36, and the error 
bar includes the distance uncertainty. The CO images show that the parent molecular cloud of 
KR 140 has been greatly disrupted by the nebula and the exciting star VES 735 and so this mass 



-15- 



will be an underestimate of the cloud's initial mass. 



7.2.1. Estimating the Mass of the Original Molecular Cloud 

One way to estimate the mass of the original molecular cloud would be to assume that 
KR 140 is indeed a blister Hn region as our data suggests. Yorke et al. (1989) have modeled 
blister Hn regions with one exciting O star and found that mass loss rate of material through the 
blister is 3 — 5 x 1CP 3 Mq yr _1 . However, these authors modeled an 06 star in a cloud with 
nn 2 = 500 cm -3 , which is not an accurate description of either VES 735 or its parental cloud. 
Unfortunately, Yorke et al. (1989) give no indication of how to scale their mass loss rate for 
different values of Q(H°) or njj 2 . A analytic estimate of the amount of mass lost from a blister Hn 
region is given by Whitworth (1979). His Equation (41) has the desirable property that it agrees 
with the results of Yorke et al. (1989) for the values used in their models; therefore, it might have 
the correct scaling for Q(H°) and n# 2 . However, his calculations were based on a cylindrical 
geometry which is not a very realistic model for a more bowl-shaped spherically symmetric region 
such as KR 140. 

In § [T^, we made use of the blister evolution formula given by Franco et al. (1994). These 
authors considered a simple symmetric blister region and were also able to estimate the cloud 
evaporation rate: 

M « 7rR 2 s m p 2n H2 C H (l + , (10) 



2Rs J 

where Rs is the Stromgren radius and m p is the proton mass. In § |5.3| , we found that an initial 
cloud density of 280 cm -3 and Rs = 0.64 pc would match the observed radius of KR 140 in a time 
of 1 Myr (consistent with the age of VES 735). Using these values in the above equation gives a 
mass loss rate of 2.9 x 10 -4 Mq yr _1 . Therefore, about 290 Mq of molecular material has been 
eroded from the cloud after 1 Myr. Including the effects of dust lowers Rs, increases the age, and 
lowers M, resulting in a very similar amount of erosion (280 Mq). 

Taking into account the mass lost through erosion and the present ionized mass, we estimate 
the mass of the original molecular cloud in which VES 735 formed to be about 4.9 x 10 3 Mq, 
so this cloud would be classified as a dwarf molecular cloud (Elmegreen 1985). This seems to 
be a remarkably small cloud mass from which an O star has formed. Williams & McKee (1997) 
performed a statistical analysis of many nearby OB associations and molecular clouds, and 
found that it was more likely that clouds with masses greater than 10 5 Mq should form O stars 
than clouds with lower masses. Other Hn regions have been found with small molecular cloud 
masses (Hunter et al. 1990), but O star formation with a cloud mass below 10 4 is considered rare 
(Elmegreen 1985). Of course we have been attracted to this cloud because of the effect of VES 735 
and not to other parts of the W3 cloud complex (Heyer et al. 1998) which have not formed O 
stars. 
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7.3. Interpretation of the Velocity 

Stepping through the velocity channels in the CO data cube near the velocity of KR 140 
shows that at higher (least negative) velocities (~ —40 km s" 1 ) there is no emission, then the CO 
emission comes in at the southern end of the Hn region and spreads northward, before filling in 
the 1420 MHz contours at a velocity of —48.78 km s _1 . Perhaps the best way to illustrate these 
data is by examining the cube in velocity-latitude space. Figure |ll] shows such a CO image of 
KR 140, averaged over the longitude range 133.366° to 133.477°. Note the "hole" at the velocity 
of the ionized gas and VES 735. 

The question arises as to the relative radial position of VES 735 and KR 140 with respect 
to the molecular gas. In general, it is difficult to use only the CO radial velocity information to 
determine absolute distances to clouds or GMCs. However, in this direction, Galactic rotation 
causes radial velocity to become more negative with increasing distance. On the face of it, that 
would place the molecular cloud on the far side of the star and nebula. However, this is possibly 
too naive an interpretation because if the cloud is only 10s of pc in radial extent (like its dimension 
in the plane of the sky), then the average shear would be too small to explain the large range in 
velocity (unless there were a large enhancement from a density wave). 

To address the question of relative position, we made use of both the Hi and CO data cubes 
to make an estimate of Ay produced by gas with velocities out to —45 km s _1 . As in § |7.2| 
we used the conversion factor X = (1.9 db 0.5) x 10 20 cm" 2 (K km s" 1 )" 1 to convert I(CO) to 
N(H2). We integrated both data cubes to obtain the atomic column densities. At the position of 
VES 735 we found N(Hi) = 4.5 x 10 21 cm -2 assuming the emission is optically thin, and N(H2) 
= 3.1 x 10 21 cm -2 , the latter mostly from local rather than Perseus arm gas. Summing the 
contributions of atomic and molecular hydrogen we obtain the total hydrogen column density, N# 
= 2N(H2) + N(Hi)= 1.1 x 10 22 cm -2 . The total visual extinction is computed using the standard 
conversion factor, Ay = 5.3 x 10" 22 mag cm" 2 (Bohlin, Savage & Drake 1978). At the position 
of VES 735 we obtain Ay = 5.7 ± 0.9, where most of the uncertainty comes from the uncertainty 
in the X factor (using the X factors of Digel et al. 1996 gives a lower Ay). This value of Ay 
compares favorably with the values around 5.7 derived by Kerton et al. (1999) using a number 
of methods (e.g., B — V colour, Ha emission measure, DIBs). However, if the cloud's molecular 
column density (summed over -45 to -53 km s" 1 ) of N(H 2 ) = 2.2 x 10 21 cm" 2 (§ |7j) is included 
in the extinction calculation, then the Ay rises to 8.0, which is quite inconsistent with the 
observed Ay to VES 735. We also made maps of the predicted extinction over the surface of the 
nebula for comparison with the extinction map derived by comparing Ha surface brightness with 
radio emission (Kerton et al. 1999). Again, inclusion of the extinction from gas in the molecular 
cloud produces more than two magnitudes too much extinction. This is strong evidence that both 
VES 735 and the ionized gas of KR 140 lie on the near side of the molecular cloud gas which has 
velocities —45 to —53 km s . 



The extinction measurements combined with Figure 



11 



suggest that KR 140 could be a 
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blister Hn region on the near side of the molecular cloud. This simple geometrical interpretation 
of KR 140 runs into difficulty if a systematic champagne flow has developed with gas flowing away 
from the parent cloud at up to the sound speed of ~ 10 km s _1 . At face value, our Ha velocity 
implies that the ionized gas is redshifted with respect to the molecular gas, which means the Hn 
region should be on the far side of the molecular cloud. However, it can be noted that Ha is 
considered to be a poor line for estimating the velocity field of champagne flows (Israel 1978; Yorke 
et al. 1984). Modeling of line profiles in champagne flows also shows that the velocity field can be 
quite random and of low amplitude, depending on geometry (Yorke et al. 1984). Recall that the 
geometry here is certainly not plane-parallel and probably more like a broken shell. Furthermore, 
our measurement sampled only gas within a projected distance of 30" of VES 735. Fabry-Perot 
data of the whole nebula in a line other than a hydrogen line would be a good way to obtain a 
better picture of the velocity structure of the ionized gas in KR 140. 

The question then arises as to the origin of the radial velocity spread within the molecular 
cloud. If this is a single, gravitationally bound cloud, then this is just the virial velocity, and its 
mass can be estimated from 

M vir = -, (11) 

where R is the length scale for the cloud, a is the linewidth, and a v i r is the virial parameter and 
includes the effects due to surface pressure, magnetic fields, and nonuniform densities (Bertoldi & 



McKee 1992). As in § 7.2 , we take R = 7 pc, and measure a typical linewidth to be 5 km s (see 



Fig. [□]). We take a vir to be 1.1 (Williams & McKee 1997). With these values M vir = 1.9 x 10 5 M Q 
which is well over an order of magnitude larger than our measured CO mass of 4.4 x 10 3 Mq. 

This result implies that the velocity field has another origin and can significantly alter the 
cloud structure over time. In fact, recalling that a velocity of 1 km s _1 corresponds to 1 pc in 
a million years, the "crossing time", 2R/a, is only 3 Myr, comparable to our above estimates of 
the age of VES 735 and the Hn region. From the radial arrangement derived from extinction, 
it seems that the background molecular material is moving toward KR 140. Perhaps the parent 
cloud giving rise to VES 735 was the result of a converging flow or "cloud collision" that is still 
ongoing (relative motion in the plane of the sky would be expected too). 



8. Hi SIGNATURE 

In theory the presence of Hll regions and their surrounding molecular material should be 
easily "visible" in Hi data sets: in both the ionized and molecular regions there should be a 
deficit of Hi emission due to a lack of neutral atomic hydrogen. In addition, if the geometry is 
favorable, neutral atomic hydrogen associated with the photodissociation region (PDR) at the 
ionized-molecular interface should be visible. In practice it is actually very difficult to observe 
Hi features that are unambiguously associated with an Hn region in raw Hi data sets. This is 
primarily due to velocity confusion along the line of sight caused by the "turbulent" motion of the 
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atomic gas. Since the velocity of the neutral atomic gas can be many times the typical channel 
width in the Hi data sets velocity becomes a much poorer proxy for physical distance than in a 
CO data set, which is probing a species with lower velocity dispersion (heavier, cooler, and less 
turbulent). The situation is especially problematic for Hn regions in the Galactic plane where 
one has to look through a large column of atomic gas towards the Hn region. For KR 140 we are 
looking through the local arm and part way into the Perseus arm of the Galaxy. A preliminary 
reconnaissance of the Hi cube confirmed that the complex Hi emission structure along the line of 
sight makes seeing any Hi signal associated directly with KR 140 extremely difficult. Nevertheless, 
some simple processing of the Hi data cube does bring out some features associated with the 
region. 

In order to exclude local Hi emission (which is assumed to have a relatively smooth spatial 
structure) and to enhance the dynamic range of the resulting channel maps we constructed a 
median-subtracted data cube (Joncas et al. 1992; Joncas et al. 1985). In this technique a median 
spectrum is calculated for the data cube and then subtracted from each spectrum making up the 
cube. The resulting channel maps thus can contain both negative and positive values indicating 
deviations relative to the median base level. 



Figure 12 shows channel maps of the median-subtracted cube over the velocity range —43.40 
to —54.95 km -1 bracketing the ionized gas velocity of —46 ± 2.0 km s _1 (Kerton et al. 1999). We 
could not detect any features definitely associated with KR 140 in the channel maps outside of 
this velocity range. 

Examining these maps we note the following three features. First, there is a noticeable deficit 
of Hi seen in velocity channels —46.70 to —51.65 km s _1 outside of the Hn region. There is excellent 
positional agreement between this deficit in Hi and the observed position of the CO emission. The 
deficit can be simply interpreted as being caused by a lack of Hi emission in the molecular material 
surrounding KR 140. Second, we also see a drop in the Hi emission occurring within the Hn 
region in the —43.40 and —45.05 km s _1 channels. This deficit is most likely associated with the 
ionized gas in KR 140 as suggested by their spatial correspondence. Figure [D| presents spatially 
averaged Hi and CO spectra for an area just outside of the Hn region to the north-east and the 
area inside the Hn region. The anticorrelation between Hi emission and the presence of CO and 
Hn is evident. Third, there is an enhancement of Hi emission at (133. 36°, 0.1583°), seen best in the 
velocity channel —46.70 km s . This could be low velocity dispersion material associated with 
the PDR. Atomic material in the PDR is expected to have a low velocity and thus should be seen 
in channels corresponding to CO emission. 



9. THE ENVIRONMENT OF KR 140 
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9.1. Spontaneous Massive Star Formation 

In the context of the Perseus Arm star formation activity, the KR 140 complex seems to be 
unique. Figure [l] shows that KR 140 is isolated from the massive and violent star formation that 
is ongoing around it. This isolation is evidence to us that KR 140 is an example of massive star 
formation in our Galaxy that is untriggered, at least in the sense used in the context of sequential 
star formation. In none of our data sets does there appear any evidence for a trigger of the star 
formation in KR 140. The exciting cluster of W4, OC1 352, is about 60 pc away from KR 140 for 
a cluster distance of 2.35 kpc (Massey, Johnson & DeGioia-Eastwood 1995). For a sound speed 
of 0.6 km s _1 (isothermal speed in H2 at 100 K), the time for a signal to reach KR 140 would 
be about 90 Myr, much greater than our estimated age of KR 140 of a few million years or the 
age of OC1 352. From W3, a signal would take about 60 Myr, but W3 is itself much less than 
10 5 years old (Kawamura & Masson 1998). The supernova remnant HB3 is also in this complex 
(Normandeau et al. 1997), but from Figure [l], the edge of the remnant is nowhere near to KR 140. 
Thus, KR 140 does not seem to be triggered by a neighboring H11 region or a nearby supernova 
remnant, unless the impulse came along the line of sight. This conclusion is complemented by the 
overall smoothness of the KR 140 nebula; while the nebula does show density inhomogeneities, it 
is not far removed from a circular shape. Thus, any perturbation that might have triggered the 
star formation within KR 140 must have been a large scale phenomenon with a characteristic scale 
of ~ 10 pc. This kind of triggering might be more consistent with triggering via a spiral density 
wave (e.g., Elmegreen 1994, 1995) or by colliding molecular clouds. Of course, we cannot rule out 
those kinds of triggers, but the observational evidence would be difficult to find. 

It is interesting to contrast KR 140 to another star forming region that has been studied 
with multiwavelength data, the Gemini OBI molecular complex (Carpenter, Snell & Schloerb 
1995a, 1995b). Within the molecular complex, these authors find young star clusters (from near 
infrared data) and a number of dense cores (as identified by CS observations) associated with 
IRAS point sources. Carpenter et al. suggest that the arc-shaped morphologies of these cores have 
been formed by swept-up gas from expanding H11 regions, and that they would form the massive 
star clusters in the region. The other lower infrared luminosity sources (most likely belonging to 
lower mass cores) in the Gem OBI complex are not found to be correlated with any arc-shaped 
structures or filaments in the molecular gas, and they are not adjacent to any H11 regions; in fact 
they are spread almost randomly around the complex. Carpenter et al. therefore conclude that 
induced star formation is the prominent mode of formation for massive stars in the Gem OBI 
molecular complex. 

If KR 140 is indeed untriggered, then it seems to be an unusual form of spontaneous star 
formation since current observations suggest that the isolated mode of star formation is generally 
associated with low mass star forming regions as seen in Gem OBI or even in Taurus. 
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9.2. IRAS Point Sources and Protostars 

As has been seen in other sites in the Galaxy, an Hn region can initiate star formation via 
the expansion of its ionization front (e.g., Elmegreen & Lada 1977). Other stars might also have 
formed spontaneously. Evidence for other star forming regions within KR 140 may be sought by 
examining the IRAS point source catalog. In addition to the IRAS point source 02165+6053 that 
is cross-referenced to VES 735, there are six other IRAS point sources in the area in and around 



KR 140 (see Table y). The crosses in Figure 14 show the positions of these six IRAS point sources. 
The circled crosses are the sources discussed here that seem most likely to be associated with the 
KR 140 complex. 

The source IRAS 02160+6057 identified with the north-west dust arc has been the subject 
of two molecular line investigations. Wouterloot & Brand (1989) identified it as a potential 
star- forming area via its IRAS colors (see their paper for the exact selection criteria), and 
examined it (along with about 1300 other IRAS point sources) for CO emission. They found a 
CO feature in that direction at a velocity of —49.7 km s _1 (LSR), which corresponds to CO in 
the associated background molecular cloud, and assigned it the catalog number [WB89]417. This 
point source was then observed in an H2O line by Wouterloot, Brand Sz Fiegle (1993), but they 
were unable to detect any emission. This line of sight has one one of the highest column densities 
in the KR 140 Hn region, and so it is possible that a protostar could be forming there as a result 
of the expansion of the Hn region. However, examination of the HIRES images shows no point-like 
features in the dust arc, and none are found in follow-up submillimetre observations with SCUBA 
(Kerton et al. 2000). Therefore, IRAS 02160+6057 is most likely simply part of the KR 140 dust 
shell. Likewise, IRAS 02168+6052 appears to just be part of the eastern side of the dust shell; it 
has not been the subject of any molecular line observations. 

The source IRAS 02157+6053 also seems to have been an identification of some part of the 
dust shell. In the submillimetre there is more structure in this region; this object appears to be a 
molecular core rather than a protostar (Kerton et al. 2000). 

The point source IRAS 02171+6058 is identified with the dust feature to the north of the 
KR 140 complex. It was included in a CS(2-1) survey by Bronfman, Nyman & May (1996) of 
IRAS point sources that have colors characteristic of ultracompact Hn regions; however, they 
were unable to make a detection. Lyder & Gait (1997) observed this source along with other 
ultracompact Hn region candidates in a search for methanol maser emission. Again, they were 
unable to detect any maser emission from IRAS 02171+6058. These non-detections do not rule 
out the possibility that the source is a protostar. In fact there is a submillimetre detection and all 
evidence seems consistent with a B5V star (Kerton et al. 2000). It is outside the radio contours of 
Hn region and is therefore difficult to interpret as being triggered. 

A photographic survey of the W3 and W4 region turned up Bright InfraRed Stars (BIRS; 
Elmegreen 1980) which are brighter in I band than in R band. There are five of these stars in the 
KR 140 region of the sky. Table gives their coordinates along with their R and I magnitudes and 
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the positions are indicated by triangles on Figure 14. Elmegreen (1980) estimated that if these 
stars were at the distance of the Persues arm then they might be deeply embedded massive early 
or pre- main sequence stars, or even giants or supergiants. From examining the overall distribution 
of the BIRS stars (Figure 1 in Elmegreen 1980), we find there is an overabundance of BIRS around 
KR 140 as compared to the average. It is difficult to say whether or not these BIRS stars are 
physically associated with the KR 140 complex. However, BIRS 128 on the eastern side of the 
nebula seems to lie at a well defined edge of the dust shell, and within the CO shell. 



9.3. The Gas Cloud 

We have been able to estimate some properties of the original molecular cloud which spawned 
the KR 140 Hn region. Originally having M ~ 4.9 x 10 3 Mq this cloud would then be classified 
as a dwarf molecular cloud. These are quite numerous in the ISM. The column density of 
2.2 x 10 21 cm" 2 is about the same as the median column density of 2 x 10 21 cm~ 2 in the molecular 
gas in the outer Galaxy (Heyer et al. 1998). Once VES 735 formed in the molecular cloud, it had 
a tremendous impact on the environment around it. It has incorporated 25 Mq, ionized about 
160 Mq, and a further 290 Mq of cloud material has been removed since the ionization front 
broke out of the front part of the molecular cloud. The remnant cloud is of fairly low density, 
tih 2 ~ 100 cm -3 . 

The velocity width of the molecular cloud is too large for the cloud to be bound, but it is 
unlikely to arise from velocity shear due to differential Galactic rotation. It is possible to interpret 
the velocity structure as a large scale flow of material towards KR 140. If true, this flow could 
have had an impact on the star formation history of this region. 

We can use our estimate of the original molecular cloud mass to estimate the mass of the 
stellar cluster that likely formed along with VES 735. Hunter et al. (1990) observed clouds 
in our mass range and estimated cluster masses using a Miller-Scalo initial mass function 



(IMF; Miller & Scalo 1979). Using their data we find the correlation shown in Figure 15: 
logM duster = (0.254 ± 0.111) log M c id + (2.306 ± 0.449). With the above cloud mass we estimate a 
cluster mass of about 1.7 x 10 3 Mq. Of course, since we know the mass of the most massive star 
we can use a Miller-Scalo IMF to estimate the mass of the cluster directly. Using Equation 12 in 
Elmegreen (1983) we estimate a cluster mass of about 1.6 x 10 3 Mq, consistent with the above 
"empirical" estimate. Using these cluster and cloud masses we estimate a star formation efficiency 
(SFE = M duster /{M c i us ter + ^cid)) of about 25% which is typical for clouds of this size (Hunter et 
al. 1990), although values can range over two orders of magnitude from cloud to cloud (Williams 
& McKee 1997). This low value of the SFE implies that the cluster within KR 140 will not be 
gravitationally bound. 

Near infrared observations (such as ones in the ongoing 2MASS survey) should be able to 
detect a number of pre-main sequence stars from this cluster around VES 735. It would be 
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interesting to compare the IMF in this region to the IMF in regions where the star formation was 
triggered. Which, if either, is a good description of the IMF in the Galaxy? 

10. CONCLUSIONS 

We have utilized our multiwavelength data set (all at a resolution of about 1') to study not 
only the physics of the Hn region itself, but also, since the data are from the larger CGPS survey, 
to study KR 140 in the context of the overall picture of star formation in the Perseus spiral arm. 
We find no evidence for a mechanism that triggered the formation of the 08.5V(e) star, VES 735, 
and its (largely unseen) cluster. We therefore conclude that this region formed spontaneously out 
of its parent molecular cloud, independent of the more vigorous star formation in W3 and W4 
nearby 

Our data of KR 140 are consistent with the model of a bowl-shaped region viewed close to 
face on. Extinction measurements to the exciting star, VES 735, and nebula show that the Hn 
region is quite likely on the near side of its molecular cloud. We have not observed any champagne 
flow, and cannot rule out other geometries. KR 140 has an age less than a few million years. We 
have estimated that the original molecular cloud had a mass of 4.9 x 10 3 M and an average 
density about 100 cm -3 , which classifies it as a dwarf molecular cloud. This makes KR 140 even 
more unusual as it is a rare example of an O star that has formed in a cloud with a mass less than 
10 4 Mq. There is tentative evidence that the molecular material is undergoing a large scale flow 
towards KR 140. Follow-up observations are needed to pursue this idea. There are four IRAS 
point sources associated with the KR 140 complex, one of which is a possible protostar candidate 
and another a molecular core. Near infrared observations of KR 140 are needed to find and study 
the young cluster that likely formed along with VES 735. 
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Fig. 1. — KR 140 in context. This 1420 MHz continuum image of the CGPS pilot region survey 
shows KR 140 (circled) in relation to the W3/4/5 star forming regions and the HB3 supernova 
remnant. The greyscale is linear with a white value of 2.0 K and a black value of 15.0 K. 
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Fig. 2. — (a) 12 fim HIRES image convolved to 1420 MHz resolution, and then overlaid with 
1420 MHz contours. The contour values are the same as Figure 3a. The greyscale is linearly 
scaled from 2.0 to 20.0 MJy sr" 1 (white-black), (b) 25 (im HIRES image convolved to 1420 MHz 
resolution, and then overlaid with 1420 MHz contours. The contour values are the same as Figure 
3a and the image is scaled from 4.0 to 25.0 MJy sr _1 (white-black), (c) Original 60 fxm HIRES 
image. The beam shape is shown in the lower right-hand corner. The image is scaled from 10.0 to 
185.0 MJy sr^ 1 (white-black), (d) Original 100 /jm HIRES image. The beam shape is displayed 
in the lower right-hand corner. The image is scaled from 40.0 to 350.0 MJy sr _1 (white-black). 
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Fig. 3. — (Top) Contours of 21 cm brightness temperature of KR 140 at 5.0, 6.0, 7.0, 8.0 (dashed), 
9.0, 9.5 (heavy), and 10.0 K. The source is slightly more extended than indicated by the lowest 
contour; the Hn region fades into the background at contour level about 3 K. The beam shape is 
displayed in the lower-right corner. (Bottom) Greyscale with linear brightness scale (1.868 — 11.0 K 
white-black). The point source at I = 133.585° and b = 0.058° is probably extragalactic as it has 
a non-thermal spectral index. The position of the exciting star of this region, VES 735, is denoted 
by the star symbol. 
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Fig. 4. — \og{L/Lo) vs. \ogT e ff with evolutionary models of Schaerer &: de Koter (1997) shown 
for three high-mass stars. Individual points in the models (crosses) are labeled with the stellar age 
(10 6 years). The point with error bars represents VES 735, the exciting star of KR 140. 
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Fig. 5. — Azimuthally averaged dust temperature profile for KR 140. Origin is located at 
I = 133.43°, b = 0.068°. The extent of the ionized region is indicated by the dashed vertical 
line. The temperature distribution was sampled every 1' in radius and every 10° azimuthally. 
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Fig. 6. — Tioo map overlaid with the 1420 MHz contours. The greyscale is linear with a white value 
of 0.0 and a black value of 0.002. Both images and the contours are at the 100 /im resolution. 
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Fig. 7. — The HIRES image of KR 140 at 25/im. The greyscale is linearly spaced with a white 
value of 4.0 MJy sr^ 1 and a black value of 30.0 MJy sr -1 . The star symbol near the center of the 
nebula corresponds to the position of the ionizing star VES 735. The beam shape of this HIRES 
map is shown in the lower right-hand corner. 
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Fig. 8. — A Cloudy model (solid line) using VES 735 stellar parameters and f c = 0.44. Integrated 
IRAS fluxes for KR 140 are shown by the crosses. A single temperature (28.25 K) v 2 B u spectrum 
is shown for contrast (dashed line). 
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Fig. 9. — An image of integrated CO emission between —45.53 km s _1 and —47.16 km s _1 (LSR). 
1420 MHz contours are overlaid with values that are the same as Figure || The image is linearly 
scaled from 0.0 to 15.0 K (white-black). 




Fig. 10. — Integrated CO emission (—45.53 to —47.16 (LSR)) associated with KR 140 with contours 
of HIRES 12 jum emission, both at ~ l' resolution. The CO image is linearly scaled from 0.5 to 
10.0 K (white - black), the contours are spaced 2 MJy/sr apart from 7.0 to 25.0. 
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Fig. 11.— Velocity-latitude CO image of KR 140 averaged over the range 133.366° < I < 133.477°. 
The image is linearly scaled from 0.0 to 7.0 K (white-black). The rectangle denotes the area where 
the ionized gas is thought to reside. The latitude range was derived from the radio continuum 
image, and the velocity range is derived from the measured radial velocity of the ionized gas 
(-46 ±2.1 km s^ 1 ). 
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Fig. 12. — Hi channel maps of the region around KR 140 from —43.40 km s to —54.95 km s . 
The Hi maps have been median subtracted and are linearly scaled from —2 to 49 K (black - white). 
The maps have also been smoothed to 2' resolution for clarity. To show the extent of the ionized 
region we have overlaid the 5 K contour from the 1420 MHz data. 
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Fig. 13. — Spatially averaged CO (dashed line) and Hi (solid line) spectra for a region containing 
the CO emission to the northeast of KR 140 (top), and an area containing the entire Hn region 
(bottom). The CO spectra have been multiplied by a factor of two. The Hi spectra have been 
median subtracted. 
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Fig. 14. — IRAS point sources (crosses) and BIRS positions (triangles) near KR 140. Circled 
crosses indicate point sources probably associated with KR 140. The linear greyscale for the 60 
image is — 185 MJy sr" 1 (white - black). 
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Fig. 15. — The correlation between the mass of a molecular cloud and the mass of the stellar cluster 
that has formed within the cloud. Derived from data given by Hunter et al. (1990). The rms scatter 
in the relation is 0.245 dex. 
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Table 1. KR 140 Physical Properties 





Property 


Value 


General 


Distance 3, 


2.3±0.3 kpc 




Diameter 


5.7 pc 


Exciting Star a 


Name 


VES 735 




Apparent V Magnitude 


12.88 ±0.01 




Apparent B — V 


1.53 ±0.02 




Spectral Type 


08.5 V(e) 




log[Q(# )] 


- 48.5 






~ 10 5 Lq 


Radio Cont. (1420 MHz) 


Flux Density 


2.35±0.05 Jy 




Avg. Emission Measure 


2000 cm" 6 pc 




n e (rms) 


30 cm- 3 




Mmi 


160 M Q 




iog[g(ff°)] 


48.05±0.1 


Infrared (IRAS) b 




110 Jy 




F 25 


140 Jy 




^60 


970 Jy 




Fwo 


2300 Jy 




ibol 


1O 45 L 


Molecular Line (CO J=l-0) 


n H . 2 


~ 100 cm" 3 




Current M totaJ 


4400 M Q 



a Data from Kerton et al. (1999) 

b Fluxes measured from HIRES images of KR 140 
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Table 2. Radio Flux Factors 



Model 


Abundance 8, 


Grains b 


Temperature 


vF v (ergs s 1 cm 1 ) 


1 


Hn 


Orion 


variable 


4.2 


2 


Hn 


ISM 


variable 


3.6 


3 


Hn 


none 


variable 


4.6 


4 


no He 


Orion 


variable 


3.9 


5 


no He 


ISM 


variable 


3.3 


6 


Hn 


Orion 


fixed (7500K) 


3.3 



a Hn - average Hn abundance as reported by Ferland (1996) 

b ISM - typical high R grains; Orion - low R grains; see Ferland (1996) 
for details 
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Table 3. IRAS Point Sources Near KR 140 



IRAS 


r 


b° 


Notes 


02168+6052 


133.49 


0.07 


part of dust shell 


02160+6057 


133.37 


0.12 


part of dust shell 


02157+6053 


133.35 


0.04 


part of dust shel; molecular core 


02171+6058 


133.49 


0.17 


no CS or CH3OH detection; colours of UCH11 region 


02174+6052 


133.56 


0.09 


not associated (?) 


02156+6045 


133.39 


-0.08 


not associated (?) 
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Table 4. Bright Infrared Stars (BIRS) Near KR 140 



Star 


1° 


b° 


R 


I 


BIRS 128 


133.50472 


0.07443 


20.5 


13.8 


BIRS 129 


133.42336 


0.11588 


19.6 


13.1 


BIRS 130 


133.42205 


0.09973 


20.5 


13.8 


BIRS 131 


133.40169 


-0.04194 


18.7 


13.6 


BIRS 132 


133.40358 


-0.04767 


18.0 


13.7 



